BABAR-CONF-02 /024 
SLAC-PUB-9319 
July 2002 



Search for the Exclusive Radiative Decays B — ^ p7 and — ^ 0^7 



The BABAR Collaboration 
July 24, 2002 



Abstract 

A search for the exclusive radiative decays B — > p{770)^ and — > a;(782)7 is performed on a 
sample of 84 million BB events collected by the BABAR detector at the PEP — II asymmetric 
e~^e^ collider. No significant signal is seen in any of the channels. We set preliminary upper limits 
of i3[B0 ^ PO7] < 1.4 X 10-6, B[B+ /j+7] < 2.3 x 10"^ and ^[5° ^ ^7] < 1.2 x lO"*' at 90% 
Confidence Level. Combining these into a single limit on the generic process B p7, we find the 
preliminary limit B[B — > p'y] < 1.9 x 10"^ corresponding, to a limit of B[B pj]/B[B K*"f] < 
0.047 at 90% Confidence Level. 
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1 Introduction 



The effective flavor-changing neutral current processes B ^ pj and uj^ probe physics at 

high mass scales both within the Standard Model and within the context of possible new physics 
scenarios through the underlying b ^ d'y "penguin" transition [Q. The decays are analogous to the 
B K*j process mediated by the b ^ s'j transition. The expected rate ol b ^ dj transitions is 
suppressed by the ratio of CKM matrix elements iVj^/Visp relative to 6 — > 57 transitions. There 
has been considerable interest recently in these exclusive channels, resulting in several calculations 
of the branching fractions expected in the Standard Model, which indicate a range B[B^ P~^7] = 
(0.9 — 1.5) X 10~^ 1^. Though the theoretical uncertainties for the branching fractions remain 
large, the possibility of extracting the ratio of CKM elements iVj^/Vj^p through the ratio B[B — > 
{p/u!)y]/B[B K*y] with less uncertainty has been explored ||3|[^. The observation of B ^ pj 
and B^ — > ivy would constitute the first evidence of the b ^ dy radiative transition and is of 
considerable interest as the first step towards extracting iVtd/Vf^,! from measurements of these 
channels. Previous searches have found no evidence for these decays [||. 

2 The BaBAR Detector and Dataset 

The decay S — > p7 is reconstructed in the modes B^ — > p^y with p^ tt^tt^ and B^ with 
p'^ — > TT+vr'^ (charge-conjugate modes are implied throughout), while B^ —>■ toy is reconstructed with 
uj — > TT~^7T~7r^. The analysis uses a sample of 84 million BB events in 78fb~^ of data collected by 
the BABAR detector Q at the PEP — II collider ||^] on the T(4S') resonance ( "on-resonance" ) , and 
9.6 fb~^ of data taken 40MeV/c^ below the T(4S') resonance ("off-resonance"). The reconstruction 
uses quantities both in the laboratory and T(4S') center-of-mass (CMS) frames, where the latter are 
denoted by an asterisk. The detector response to the signal and background processes are studied 
with a detailed Monte Carlo simulation based on Geant4 Q and cross-checked with control samples 
in the data. The off-resonance data provide a control sample of the primary backgrounds from the 
continuum production of u, d, s, and c quark-antiquark pairs, while exclusively reconstructed 
B — > Dtt decays provide a sample to cross-check the simulation of BB events. 

3 Analysis Method 

3.1 Reconstruction of the Primary Photon 

The primary photon in the decay is identified as a local maximum within a contiguous deposition of 
energy in the crystal array of the electromagnetic calorimeter (EMC). We require that the photon 
lies in the calorimeter acceptance of —0.74 < cosO < 0.93, where 6 is the polar angle to the 
detector axis. The energy of the photon, measured in the center of mass system, must satisfy 
1.5 < E* < 3.5GeV. The photon candidate is required to be isolated from all other local maxima 
in the calorimeter by 25 cm. It must also be inconsistent with the trajectories of all reconstructed 
charged tracks. Photons consistent with vr*^ and r] production are vetoed by removing candidates 
that form an invariant mass within 20(40) MeV/c'^ of the T^^{r]) mass when paired with any other 
photon in the event with energy greater than 50(250) MeV. Energetic vr'^s which produce photons 
that cannot be resolved as separate local maxima are suppressed by requiring the lateral profile of 
the energy deposition to be consistent with a single photon. 
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Figure 1: Pion selection performance as measured in data using D*^ 
as a function of momentum in the laboratory frame piab- 



D°Tr+,D^ K-TI+ decays 



3.2 Reconstruction of Charged Tracks 

The charged tracks used in identifying the p/u meson are reconstructed in the silicon vertex detector 
(SVT) and drift chamber (DCH) and are required to have a trajectory consistent with production 
near the beam interaction point as well as a minimum of 12 hits in the drift chamber. A charged 
pion selection based on dE/dx measurements and Cherenkov photons reconstructed in the ring- 
imaging Cherenkov detector (DIRC) is used to reduce backgrounds from B — > K*^ and other 
b ^ S'j processes by vetoing charged kaons from these processes. Both the reconstructed Cherenkov 
angle and the number of Cherenkov photons observed is required to be consistent with the pion 
hypothesis. Figure ^ shows the particle identification performance achieved. 



3.3 Reconstruction of p/uj Meson 

The TT^TT^ candidates are reconstructed by calculating a common vertex for two tracks of 

opposite charge. For the vr+vr'^ and lo tt^tt^tt^ reconstruction, vr'^ candidates are identified 

as two photon candidates reconstructed in the calorimeter each of with energy greater than 50 MeV. 
The invariant mass of the pair is required to be 0.115 < M^^ < 0.150 GeV/c^. A kinematic fit with 
My^ constrained to the nominal vr'^ mass is used to improve the momentum resolution. The p^ 
candidates result from vr'^ candidates paired with an identified charged pion. The invariant mass 
Mtttt of the p candidates is required to be between 0.520 and 1.020 GeV/c^. The uj candidates 
are reconstructed from combinations of oppositely charged identified pions with a successfully 
calculated vertex and vr'' candidates with invariant mass 0.7596 < M^+^-^o < 0.8056 GeV/c^. 
The momentum of the candidate p mesons, measured in the center of mass system must satisfy 
2.3 < p* < 2.85 and the uj mesons must satisfy 2.4 < < 2.8. This cut, which has very high 
signal efficiency, is applied in order to reduce the number of events where more than one candidate 
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satisfies all of the cuts. 

3.4 Reconstruction of the B Meson 

The photon and p/u meson candidate are combined to form the B meson candidate. The kinematic 
properties of the B meson are evaluated in the CMS using the variables A£?* = E"^ — -E^gam 

the beam-energy substituted mass mss = \j E^eam ~ P's ' where E^gam energy of the beam, 

E'^ = E* + E*i^ is the reconstructed energy of the B meson candidate and p'^ is its momentum. 

For the purposes of the to-es calculation, p'^ is modified by scaling the photon energy so as to make 
E* + E*i^ — El^^^ = 0, under the assumption that the resolution of the primary photon energy 
dominates the /^E* resolution. This procedure reduces the tail in the mss resolution that results 
from the asymmetric photon energy response in the EMC. The signal events have = ""t-b and 
A£^* = up to the experimental resolution of ~ 3 MeV/ dominated by the beam-energy spread 
in the former case, and ~ 50 MeV dominated by the reconstructed photon energy resolution in the 
latter. 

We consider candidates in the region —0.3 < /S.E* < 0.3 GeV and 5.2 < uies < 5.29GcV/c^ and 
define a signal region of -0.2 < AE* < 0.1 GeV and 5.27 < niES < 5.29 GeV/ c^. Selection criteria 
have been optimized for best S'^/{S + B), where S and B are the expected signal and background 
yield assuming B[B^ p°7] = B[B'^ —>■ coj] = ^B[B~^ — *■ = 10~^ (as expected from isospin 
symmetry) without knowledge of the yield or distribution of events in the signal region. The signal 
region extends lower on the negative side of AE* due to the asymmetric photon energy response in 
the calorimeter resulting from energy leakage. For the small fraction of events (2.8 % for B^ 
signal) in which more than one B meson candidate satisfies all the cuts, the candidate with the 
smallest value of |A£^*| is selected. 
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3.5 Suppression of Background 

The continuum and initial state radiation backgrounds are suppressed by a neural network that 
combines event topology variables into one discriminating variable |^. The neural network responds 
non-linearly to the input variables and exploits correlations between the variables. The input 
variables are: 

• cos 9^, the cosine of the angle between the photon and the thrust axis of the rest of the event 
(excluding the B meson candidate). 

• cos^Ij, the cosine of the angle between the B meson momentum and the beam axis. 

• The energy flow in 10° bins centered on the photon candidate momentum in the CMS. 

• cos 6h, the cosine of the helicity angle. For B — > p7, we define Oh as the angle between the 
TT~^ momentum in the p rest frame and the p momentum in the B meson rest frame. For 
B^ LO'j, Oh is defined as the angle between the normal to the plane defined by the vr+vr^vr^ 
momenta in the oj rest frame and the u: momentum in the B meson rest frame, cos Oh should 
follow a sin^ Oh distribution for signal, while the continuum background is approximately flat. 

• i?2, the ratio of second and zeroth order Fox- Wolfram moments in the frame recoiling from 
the photon momentum. This is effective against initial-state radiation, since in that frame 
the jet structure of the hadrons is recovered. 

• The net flavor content, defined as \ ~ where A^^^ are the number of e^, /i^, 
and slow pions of each sign identified in the event. 

• |Az|, the vertex separation of the B meson candidate and the rest of the event along the 
beam axis, is used for B^ — > p^'^ and B^ — > 1^7. 

• cos Od, the cosine of the Dalitz angle of the to decay is used for B^ loj. cos Od is defined 
as the angle between the. tt" and the 7r~^ momenta in the rest frame of the 7r"'"7r~ system. We 
expect cos 6*0 to be uniformly distributed for the combinatorial background and to follow a 
sin^ Od distribution for true u! decays. 

A separate neural network is trained for each mode using the back-propagation algorithm on 
samples of Monte Carlo-simulated signal and background. The output of the neural network, 
defined such that the signal processes peak at one and the continuum background at zero, is cross- 
checked on an independent sample of Monte Carlo-simulated events and data control samples for 
both the signal and background. The neural network output for B^ p^j is shown in Figure |2|, 
where the Monte Carlo simulation of the continuum background is compared with the off-resonance 
data and the output for Monte Carlo-simulated B^ D~tt~^ decays is compared with events 
reconstructed in the data. This latter check gives us confidence that the Monte Carlo simulates 
the B ^ pj efficiency well, since most of the input variables of the neural network are not based 
on the properties of the signal decay itself, but rather on the properties of the other B meson in 
the event. 

We make a selection on the neural network output optimized for best S'^/ (S + B) for each mode. 
For B^ — > p'^^, an additional requirement of | cosOh\ < 0.6 is made to reject B^ p'^n^ events 
which have a cos^ Oh distribution, different from the expected sin^ Oh distribution from the signal 
process. 
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niEs (GeV/c2) 

Figure 3: /S.E* vs. ttt-es for a.) p^j, b.) — > p'^^ and c.) B^ — > toj candidates. 



3.6 Background Estimation 

Table |l] shows estimates for the expected background remaining in the signal region (—0.2 < AE* < 
0.1 GeV and mss > 5.27GeV/c^) after the neural network selection, using the ofF-resonance data 
for continuum and Monte Carlo simulation for BB backgrounds. There is good agreement between 
Monte Carlo estimates of continuum background and off-resonance data. The BB background is 
smaller than the continuum background. 

3.7 Signal Extraction 

After the neural network selection, the signal extraction in B ^ pj is performed with an unbinned 
extended maximum likelihood fit in the variables m-Es, AE* and AIj^t^ with signal and continuum 
background components. Since the BB backgrounds are small relative to the continuum back- 
ground, the signal extraction uses only a continuum component to describe the background. Biases 
due to BB backgrounds are estimated in the systematic studies described in Section ^. The signal 
mEs and AE* distributions are described by the Crystal Ball lineshape |jll|, with the exception 
of the rriES distribution for B^ — > p^-y, where the Gaussian distribution is used. Here, we expect 
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p^-i B+ p+-f B^ ujj 
(Events) 



Off-resonance data 


42. S 


,± 


7.0 


66.5 




8.8 


7.3 




3.3 


uds + cc + T+r^ Monte Carlo 


39.J 


,± 


3.8 


59.8 




4.6 


8.4 




1.9 


BB Monte Carlo 


6.6 ± 


1.0 


12.2 




3.6 


0.7 




0.3 


Signal Expectation 


9.9 ± 


0.2 


12.1 




0.3 


3.4 




0.1 



Table 1: Estimates of continuum background (using both off-resonance data and Monte Carlo sim- 
ulation), BB backgrounds and signal in the signal region. All estimates are scaled to a luminosity 
corresponding to 84 million BB pairs. Errors are statistical only. 



the photon energy rescaling to eliminate the tail. The background ttt-es and AE* distributions are 
described by the ARGUS threshold function |12| and second order polynomial, respectively. For 
B — > p'j, the Breit-Wigner lineshape is used for the signal M,r7r distribution, while the background 
incorporates a sum of a resonant Breit-Wigner component and a continuum component described 
by a first order polynomial. 

The parameters of the signal probability distributions are obtained from the Monte Carlo simu- 
lation and cross-checked in the data with the decays B^ — > K*^^, K*^ K~^-k~ for B^ p^^ and 
B^ i^*"'"7, K*^ — > K^TT^ for B^ P^l: which are topologically and kinematically similar. The 
parameters of the continuum background distributions are determined in the fit, with the excep- 
tion of the fraction of the resonant p irir contribution to the continuum M.^.,^ distribution, which 
is fixed to the value measured in the off-resonance data. In the B^ wy, the signal extraction 
is performed in a similar fit to the m-Es and AE* distributions; the M^+^-^o distribution is not 
included in the fit, because the continuum background contains a large and uncertain fraction of 
true uj decays, which is difficult to model. 

By inverting the pion selection on the charged pion (and selecting kaons) in the B^ — > p'^'y 
analysis and on one of the charged pions in the B^ p^-y analysis, we obtain an orthogonal sample 
of events enhanced in the decays 5+ i^*+7, ^ K+tt^ and B^ K*^-/, K*^ ET+vr", 
respectively. The yield of i? ^ K*'y in these two samples is determined using the same fit proce- 
dure described for B — > p^, with the expected signal distributions determined from Monte Carlo 
simulation. The resulting yields are in agreement with the expectations from previous measure- 
ments of B[B K*^] ijl^, as shown in Table thus providing a cross-check on the event selection 
and signal extraction procedure up to the statistical uncertainty in the extracted yields and the 
uncertainties in the measured branching fractions. 

The AE* vs. rn-Es distributions of the B ^ pj and B^ lo^ candidates are shown in Figure |^ 
and the fitted yields summarized in Table |3[ The quality of the fit is determined by comparing the 
minimum —logC of the fit, where C is the overall likelihood of the fit, with values obtained from 
parameterized Monte Carlo simulation and found to be in good agreement. 



4 Systematic Studies 

The systematic uncertainties in this analysis are associated with uncertainties in the efficiency of 
the signal process reconstruction predicted by the Monte Carlo simulation and with the signal 
extraction procedure. Table ^ summarizes these uncertainties. The efficiency of the track selection 
is calculated by identifying tracks in the silicon vertex detector and evaluating the fraction that is 
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Mode 


Fitted Yield (Events) 


Expected Yield (Events) 




343.2 ±21.0 


332 ± 36 




93.1 ± 12.6 


105 ± 18 



Table 2: The fitted yields in the B — > /C*7-enhanced sample described in the text, and the 
expected yield from the measured branching fractions in p!o[ |. 



Mode 




Yield 


90% C.L. Upper Limit Yield 


Bias 


Efficiency 






(Events) 


(Events) 


(Events) 


(%) 


B^^ 


p'l 


4.8 ±5.2 


12.4 


[-0.5,0.8] 


12.3 


B+ ^ 




6.2 ±5.5 


15.4 


[-0.1,2.0] 


9.2 


B^^ 


W7 


0.1 ±2.3 


3.6 


[-0.3,0.5] 


4.6 



Table 3: The fitted yields, the ranges of observed biases from BB backgrounds and selection 
efficiencies for B^ p^j, B'^ /9+7 and B^ uj^ in the on-resonance data sample. The 
efficiencies include the branching fractions for oj 7r+7r-7r° (88.8%), p° ^ vt+tt- (« 99%) and 
p+ ^ TT+vrO (f« 100%). 



well-reconstructed in the drift chamber. The pion identification efficiency in the DIRC is derived 
from a sample of D*~^ — > D^ir'^ , K~it~^ decays. The photon and -k^ efficiencies are measured 

by comparing the ratio of events N{t^ h^'K^)/N{T^ /i^vr'^vr'^), where = vr^, K^, to 
the previously measured branching ratios [13|. The photon isolation and vr'^/r/ veto efficiency are 
dependent on the event multiplicity and are tested by "embedding" Monte Carlo-generated photons 
into both an exclusively reconstructed B meson data sample and a generic BB Monte Carlo sample. 
The efficiencies of the neural network selection are compared between the exclusively reconstructed 
B — > Dtt events in the data and the Monte Carlo simulation of both the B Dtt and signal 
processes, and the observed variations taken as systematic uncertainties. The M^+^r-TrO selection 
for the B^ — > ujj is checked by comparing the resolution of the u mass peak in the on-resonance 
data with the Monte Carlo simulation and evaluating the variation in efficiency. 

Systematic uncertainties in the signal extraction procedure result from the modeling of the BB 
backgrounds and from uncertainties in the signal probability distribution functions used in the fit. 
The biases due to BB backgrounds are estimated by varying the distributions as well as the rates 
of the dominant backgrounds coming from 6 ^ 57 and B~^ p'^ir^ decays. The full range of 
biases obtained from these variations is taken as the allowed range, as shown in the column labeled 
"Bias" in Table ^. The uncertainties resulting from the fixed parameters describing the signal 
distributions in the fit are estimated by varying the parameters within the uncertainty obtained 
from the analogous B K*^ processes in the data used to cross-check the expectations from the 
Monte Carlo simulation. The effects of these variations on the fitted signal yield in each mode is 
calculated and the range of observed biases are taken as systematic uncertainties in the signal yield. 



5 Physics Results 

Based on the observed signal yields we determine 90% upper limits on the branching fraction by re- 
fitting the distributions with increasing fixed signal yields until the — logi2 deviates by 0.82 relative 
to the minimum value. We correct for bias from BB backgrounds by applying the smallest observed 
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Selection Criteria 


p% B+ p+-f B^ ujj 
Systematic Uncertainty 

(%) (%) (%) 


B Count 


1.1 


1.1 


1.1 


7 EfF. 


1.5 


1.5 


1.5 


TfO Eff. 


- 


5.0 


5.0 


TT^/f] Veto 


1.0 


1.0 


1.0 


7 Dist Cut 


2.0 


2.0 


2.0 


Tracking Eff. 


2.5 


1.3 


2.4 


TT Selection 


6.0 


3.0 


6.0 


M^+^-^o Selection 






2.0 


Neural network Selection 


8.0 


6.0 


14.0 


Fit Distributions 


5.0 


10.0 


5.0 


Total 


11.8% 


13.4% 


17.3% 



Table 4: Summary of systematic uncertainties for B ^ pj and B — > ujj expressed as percent 
error of the branching fraction. 



bias to the signal yield (increasing the signal yield in all cases). The yield of events in the data 
sample are calculated with the Monte Carlo-derived efficiency lowered by one standard deviation in 
the systematic error. Finally, the estimated number of BB events in the sample is reduced by one 
standard deviation. The resulting preliminary 90% confidence level upper limits for the branching 
fractions are B[B'^ p^j] < 1.4 x 10"^, B[B+ p+j] < 2.3 x 10"^ and ^[5° ujj] < 1.2 x 10~*^. 

For the purpose of comparing the limits to B[B K*j], we combine these limits into a single 
limit on the generic process B ^ pj defined as 

B[B pj] = B[B^ p+7] = 2 X B[B^ p^-f] = 2 x B[B° ^7]. 

as expected from isospin symmetry. The resulting preliminary 90% confidence level upper limit is 
B[B —>■ p7] < 1.9 X 10"^. Using the measured value of B[B K*'y] |^], this corresponds to a ratio 
of 

B[B pi]/B[B K*-i] < 0.047. 
We can convert this ratio to a limit on |T4(i/Vis| using the following expression from Q. 

/ I -ml/Ml V . 
\l-ml,/Mlj 

Taking the conservative limits of the two parameters, C = 0.7 and Aii = —0.25, we find 
\Vtd/Vts\ < 0.36 at 90% Confidence Level. 

6 Summary 

In conclusion, we have found no evidence for the exclusive 6 — > ^7 transitions B ^ pj and B^ — > lu^ 
in 84 million BB decays studied with the BABAR detector. The preliminary 90% confidence level 
upper limits on the branching fractions are significantly improved and within a factor of two of the 
largest Standard Model predictions, which indicate a range B[B~^ —>■ p'^'j] = (0.9 — 1.5) x 10~^. 



B[B ^ P7J 
B\B K*-y] 



Vtd 

Vts 
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